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摘 要 : 干旱 是 致 灾 最 为 严重 的 极端 气候 事件 之 一 ,研究 未 来 气候 变 暧 背景 下 旱灾 风险 的 变化 有 利 


于 科学 推进 防 灾 减 灾 工 作 部 署 。 利 用 第 六 次 国际 耦合 模式 比较 计划 的 20 个 气候 模式 数据 计算 了 
标准 化 降水 燕 散 指数 ,提取 了 基准 期 及 全 球 2%C .3 CA C 温 升 情景 下 中 国 的 干旱 特征 变量 并 计算 


干旱 危险 性 指数 ,基于 承 灾 体 预 估 数 据 计算 干旱 暴露 度 指数 和 干旱 脆弱 性 指数 ,综合 计算 旱灾 风 
险 指数 ,从 而 分 析 中 国 旱 灾 风 险 分 布 格局 并 基于 地 理 探 测 器 对 未 来 旱灾 风险 变化 进行 空间 归 因 分 
析 。 结 果 表明 :干旱 危险 性 指数 .干旱 骏 露 度 指 数 和 干旱 脆弱 性 指数 的 空间 分 布 分 别 表现 为 西北 


和 东南 相对 较 高 、 东 高 西 低 、 


西高 东 低 ;旱灾 风险 指数 具有 东 高 西 低 的 分 布 特点 ,呈现 以 高 值 集聚 
和 低 值 集聚 为 主 的 空间 正 相 关 ; 随 着 温 升水 平 的 升 高 ,未 来 旱灾 风险 以 增加 为 主 ,东部 沿海 地 区 增 


加 最 为 明显 ;人 口 数量 变化 .CDP 变 化 和 耕地 点 比 变化 是 影响 旱灾 风险 变化 的 主导 因素 。 
X 键 词 : 旱灾 风险 ; CMIP6; 温 升 情景 ; 空间 自 相关 ; 地 理 探测 器 
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在 全 球 变 暖 背景 下 ,近年 来 世界 各 地 的 极端 天 气 
气候 事件 (高 温 热 浪 .干旱 、 强 降水 等 ) 日 益 频 繁 "。 
联合 国政 府 间 气候 变化 专门 委员 会 (Intergovern- 
mental Panel on Climate Change , IPCC ) 第 六 次 评估 报 
告 (the Sixth Assessment Report, AR6) 指 出 , 2011— 
2020 年 全 球 平均 气温 较 工 业 化 前 水 平 升 高 了 约 
1.09 ,未 来 全 球 温 升 水 平 的 持续 增加 将 导致 危害 
多 发 并 发 "。 干旱 作为 一 种 频率 高 持续 久 、 人 危害 大 
的 极端 气候 事件 ,对 水 资源 供应 .农业 生产 .生态 保 
护 等 诸多 方面 产生 负面 影响 。 中 国 是 一 个 干旱 灾 
害 频 发 的 农业 大 国 ,2017 一 2021 年 全 国 年 均 作物 受 
量 面积 为 7.8x10' hm ,旱灾 造成 的 年 均 直 接 经 济 损 
失 达 3.1x10" 元 ”。 研 究 表明 ,受气 候 变 暖 影响 ,过 
去 几 十 年 中 国 整体 呈 干 旱 化 趋势 ,21 世 纪 中 后 期 变 
干 趋势 将 持续 ”, 危 及 经 济 社会 可 持续 发 展 。 因 此 ， 
科学 评估 未 来 旱灾 风险 显得 极为 重要 ,能 够 为 旱灾 
风险 预警 和 防范 提供 科学 文 返 。 

旱灾 风险 指 干旱 的 发 生 对 社会 经 济 和 生态 环 
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境 造 成 不 利 影响 的 可 能 性 ,是 致 灾 因子 危险 性 、 承 
灾 体 暴露 度 和 脆弱 性 相互 作用 的 结果 <。 旱灾 的 致 
灾 因 子 是 指 气象 系统 异常 造成 的 干旱 ,干旱 危险 性 
分 析 一 般 是 基于 干旱 指数 识别 干旱 过 程 ,提取 用 于 
定量 描述 干旱 事件 的 干旱 特征 变量 (干旱 历时 、 干 
旱 强度 等 ) 并 加 以 分 析 “。 旱 灾 的 承 灾 体 为 受 旱 灾 
影响 和 威胁 的 人 类 社会 主体 ,在 同等 致 灾 因 子 作用 
下 , 承 灾 体 的 暴露 度 和 脆弱 性 越 大 ,风险 越 大 。 暴 
露 度 反 映 了 暴露 在 干旱 下 的 人 口 和 资产 的 规模 ,可 
用 人 口 密 度 .CDP 等 指标 进行 量化 ;脆弱 性 反映 了 
承 灾 体 对 干旱 的 敏感 性 和 应 对 能 力 , 可 用 人 均 
GDP 土地 利用 情况 等 指标 进行 量化 。 不 少 研究 从 
干旱 危险 性 ` 承 灾 体 暴露 度 和 脆弱 性 3 个 方面 进行 
分 析 ,构建 3 个 综合 指标 并 计算 其 乘积 以 表征 旱灾 
风险 ,为 旱灾 风险 评估 框架 的 构建 提供 了 参考 。 
目前 ,旱灾 风险 评估 多 基于 历史 资料 开展 ,而 
通过 引入 未 来 情景 数据 可 开展 未 来 旱灾 风险 预 
估 。 世 界 气候 研究 计划 (WCRP) 组 织 的 国际 耦合 模 
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式 比 较 计 划 已 进入 第 6 阶段 (CMIP6) ,所 提供 的 气 
候 模 式 数 据 在 未 来 气候 预测 研究 中 应 用 广泛 一 。 
同时 ,基于 共享 社会 经 济 路 径 (SSPs ) 的 社会 经 济 要 
素 ( 人 口 经 济 .土地 利用 等 ) 定 量 模拟 和 预 估 数 据 
集 逐 渐 丰 富 ,为 气象 灾害 的 承 灾 体 暴露 度 和 脆弱 性 
评估 提供 了 数据 文 撑 ”。 已 有 研究 结合 气候 模式 
数据 和 承 灾 体 预 估 数 据 对 全 球 和 区 域 尺度 的 未 来 
旱灾 风险 进行 了 预 估 " 六 ,但 对 不 同 温 升 情景 关注 
较 少 。 为 了 更 好 地 应 对 全 球 温 升 加 剧 对 旱灾 应 急 
管理 带 来 的 挑战 ,有 必要 预 估 未 来 不 同 温 升 情景 
干旱 对 经 济 社会 的 影响 。 综 上 ,本 研究 基于 CMIP6 
气候 模式 数据 和 承 灾 体 预 估 数 据 评估 2% .3 % .4%C 
温 升 情景 下 的 中 国 旱 灾 风 险 及 其 变化 ,以 期 为 全 球 
变 暖 背景 下 抗旱 减灾 方案 制定 提供 科学 依据 。 


1 数据 与 方法 


1.1 研究 区 概况 

中 国 地 处 亚洲 东部 ,是 一 个 地 域 广阔 、 地 形 复 
ZR .气候 多 样 的 国家 。 为 对 比分 析 不 同 区 域 的 旱灾 
风险 特征 ,本 文 根 据 前 人 研究 将 中 国 划分 为 7 个 自 
然 地 区 , 即 东北 湿润 羊 湿润 温带 地 区 、 华 北 湿润 
半 湿 润 暧 温带 地 区 、 华 中 华南 湿润 亚热带 地 区 、 华 
南 湿润 热带 地 区 .内 蒙 草原 地 区 PE AE Aie YR He XAN 
青藏 高 原 地 区 ,分 别 用 A、B、C.D、EF 和 G 表 示 ( 图 
1)。 该 区 划 方 案 按照 地 理 位 置地 貌 . 气 候 、 植 被 等 


0  500km 
图 例 

一 国界 -- 未 定 国界 一 区 划 线 
自然 地 区 


EDI 华南 湿润 热带 地 区 
—A 东北 湿润 半 湿 润 温 带 地 区  —— E 内 蒙 草原 地 区 
华北 湿润 半 湿 润 暧 温带 地 区 ES PUO DERE 
TE 华中 华南 湿润 亚热带 地 区 BE PRSE 


审 图 号 : GS(2023)4315 号 


图 1 中 国 7 个 自然 地 区 


Fig. 1 Seven natural regions of China 


因素 的 差异 性 和 相似 性 进行 划分 所 ,已 广泛 应 用 于 
中 国 干旱 相关 的 研究 1。 

1.2 数据 来 源 

1.2.1 和 气象 观测 数据 CN05.1 数 据 集 基于 中 国 2400 
多 个 气象 台 站 1961 年 至 今 的 观测 资料 使 用 距 平 副 
近 法 插值 得 到 ,分 辨 率 为 0.25°, 是 日 前 中 国 区 域 最 
精确 的 格 点 化 气象 观测 数据 集 “”。 本 研究 将 
1961 一 2014 年 逐 月 气象 要 素 ( 降 水 量 .最 高 气温 最 
低 气温 平均 风速 ,相对 湿度 和 日 照 时 数 ) 双 线性 插值 
为 0.5° 分 辨 率 ,用 于 评估 CMIP6 历 史 气 候 模 拟 数据 。 
1.2.2 气候 模式 数据 CMIP6 气 候 模 式 数据 由 WCRP 
提供 (https://esgf-node.linl.gov/projects/cmip6) ,是 气候 
变化 预测 研究 中 最 先进 、 最 可 靠 的 数据 源 之 一 "1。 
历史 气候 模拟 试验 提供 了 1850 一 2014 年 的 历史 气 
候 模拟 数据 ,情景 模式 比较 计划 提供 了 2015 一 2100 
年 的 未 来 气候 预 估 数 据 ”。 本 研究 使 用 1995 一 
2100 年 逐 月 的 平均 气温 数据 计算 全 球 温 升水 平 ,使 
用 1961 一 2100 年 逐 月 的 降水 量 .最 高 气温 .最低 气 
温 、 风 速 .相对 湿度、 地 表 下 行 短波 辐射 .气压 数据 
计算 干旱 指数 。 情 景 模式 比较 计划 设计 的 气候 预 
佑 情景 是 不 同 SSPs 与 辆 射 强迫 水 平 的 组 合 情 景 ,本 
研究 选取 了 SSP5-8.5 情 景 (化 石 燃料 驱动 的 发 展 路 
径 与 高 辐射 强迫 情景 ) ”, 这 是 唯一 能 够 达到 4 "Ci 
升水 平 的 气候 预 估 情景 。 选 取 了 包含 所 需 变量 的 20 
个 气候 模式 ( 表 1) ,将 各 模式 数据 双 线性 插值 为 0.$。 
1.2.3 高 程 数 据 美国 地 质 勘探 局 提供 了 2010 年 全 
球 多 分 辨 率 地 形 高 程 数据 (GMTED2010) ^" ,可 满足 
不 同 尺 度 和 不 同 领域 的 应 用 需求 。 本 研究 将 30" 分 
辩 率 的 平均 高 程 数据 采用 均值 聚合 方法 重 采 样 为 
0.5" 分 辨 率 , 用 于 计算 干旱 指数 。 

1.2.4 人 口 和 GDP 数 据 姜 彤 等 ”发布 了 5 种 SSPs 
的 人 口 和 经 济 预 佑 数据 ,包含 2010 一 2100 年 0.5° 分 
辩 率 的 中 国 逐 年 人 口 和 GDP 数 据 。 该 数据 集 考虑 
了 人 口 政策 变化 对 人 口 结构 和 经 济 的 影响 ,比较 贴 
近 中 国 国 情 ”。 本 研究 使 用 SSP5 的 人 口 和 GDP 数 
据 进行 承 灾 体 暴露 度 和 及 弱 性 评估 。 

1.2.5 土地 利用 数据 土地 利用 协调 (LUH2) 数 据 集 
是 CMIP6 试 验 的 基础 驱动 数据 ,包含 850 一 2100 年 
0.25° 分 辨 率 的 全 球 逐 年 土地 利用 数据 ”i。 本 研究 
基于 历史 数据 和 SSP5-8.5 情景 数据 提取 1995 一 
2100 年 中 国 的 耕地 占 比 和 灌溉 占 比 , 即 每 个 格 点 的 
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表 1 研究 所 使 用 的 20 个 CMIP6 气候 模式 信息 
Tab.1 Information of 20 climate models from CMIP6 used in this study 


序号 模式 名 称 机 构 简称 国家 经 纬 向 格 点 数 
1 ACCESS-CM2 CSIRO-ARCCSS 澳大利亚 192x144 
2 ACCESS-ESMI-5 CSIRO 澳大利亚 192x145 
3 CanESM5 CCCma 加 拿 大 128x64 
4 CanESMS-1 CCCma 加 拿 128x64 
5 CAS-ESM2-0 CAS 中 国 256x128 
6 CMCC-ESM2 CMCC 意大利 288x192 
7 EC-Earth3 EC-Earth-Consortium 欧盟 512x256 
8 EC-Earth3-CC EC-Earth-Consortium 欧盟 512x256 
9 EC-Farth3-Veg EC-Earth-Consortium 欢 盟 512x256 
10 EC-Earth3-Veg-LR EC-Earth-Consortium 欧盟 320x160 
11 FGOALS-g3 CAS "mg 180x80 
12 FIO-ESM-2-0 FIO-QLNM 中 国 288x192 
13 GFDL-ESMA NOAA-GFDL 美国 288x180 
14 INM-CM4-8 INM 俄罗斯 180x120 
15 INM-CM5-0 INM 俄罗斯 180x120 
16 IPSL-CM6A-LR IPSL 法 国 144x143 
17 MIROC6 MIROC 日 本 256x128 
18 MPI-ESMI-2-HR MPI-M 德国 384x192 
19 MPI-ESMI-2-LR MPI-M 德国 192x96 

20 MRI-ESM2-0 MRI 日 本 320x160 


耕地 面积 和 汶 溉 面积 占 该 格 点 面积 的 比例 ,采用 面 
积 加 权 平 均 方 法 将 分 辩 率 统一 为 0.3" ,用 于 承 灾 体 
AE t BERG SS TER TER; 

13 研究 方法 

1.3.1 素 勒 图 方法 泰勒 图 方法 是 评估 气候 模式 模 


间 段 。 
1.3.3 旱灾 风险 评估 模型 

(1) 干旱 危险 性 分 析 

使 用 1 个 月 时 间 尺 度 的 标准 化 降水 蒸 散 指数 
(SPED 作为 监测 气象 干旱 的 指标 ,无 旱 、. 轻 旱 .中 


拟 性 能 的 一 种 常用 方法 ,将 相关 系数 (r) 标准 差 
(SD) 以 及 均 方 根 误差 (RMSE)3 个 统计 量 显示 在 一 
张 图 中 , 较 直 观 地 给 出 模拟 值 与 观测 值 的 差异 ”。 
本 人 研究 使 用 标准 化 泰勒 图 ,将 观测 值 和 模拟 值 的 SD 
和 RMSE 除 以 观测 值 的 SD, 以 消除 量 纲 。 标 准 化 泰 
WEP ,观测 值 的 SD 为 1,RMSE 为 0。 模 拟 性 能 较 
好 的 气候 模式 其 -和 SD 接近 1 有 日 RMSE 接 近 0。 

1.3.2 不 同 温 升 情景 的 确定 根据 IPCC AR6,1995 一 
2014 年 的 全 球 变 暧 观测 值 为 0.85 WC”。 因 此 ,本 研 
究 将 1995 一 2014 年 作为 温 升水 平 计 算 和 旱灾 风险 
评估 的 基准 期 。 基 于 CMIP6 各 模式 逐 月 的 平均 气 
温 格 点 数据 ,采用 面积 加 权 平 均 法 计算 1995 一 2100 
年 的 逐年 全 球 平均 气温 。 对 于 每 个 模式 ,计算 其 20 a 


te GER 、 特 旱 对 应 的 SPEI 数 值 范围 分 别 为 :>-0.5、 
(1.0,-0.5] (=-1.3,-1.0] (2.0, -1.5] ,«- 2.0, 
SPEI 考 虑 了 莱 散 量 对 气温 变化 敏感 的 特点 ,该 指 
数 由 降水 量 与 潜在 蒸 散 量 (PET) 的 差 值 构建 ,计算 
方法 可 参考 国家 标准 《气象 干旱 等 级 》”。 本 人 研究 
中 PET 的 计算 采用 Penman-Monteith 方 法 ,PET 和 SPEI 
的 计算 使 用 RR 语言 SPEI 包 ”。 首 先 基 于 1961 一 
2014 年 数据 计算 各 模式 SPEI, 再 将 所 得 拟 合 参数 用 
于 2015 一 2100 年 SPEI 的 计算 。 

参照 前 人 研究 喇 ,基于 SPEI 采 用 三 阔 值 游程 理 
论 方 法 识别 干旱 事件 ,选取 的 3 个 赋值 分 别 为 0.-0.5 
和 -1。 提 取 基 准 期 和 2% .3 % .4% 温 升 情景 (后 续 
简称 4 种 情景 ) 的 20 a 期 间 干 旱 频 次 干旱 历时 和 干 
旱 强 度 3 个 特征 变量 。 其 中 ,干旱 频次 是 平均 每 年 


滑动 平均 气温 ,与 基准 期 平均 气温 作 差 ,将 温差 
首次 达到 1.15 CC、2.15 CC、3.15 MC 的 未 来 某 个 20 a 
时 间 段 定义 为 2% 3% 、 4s% 温 升 情景 对 应 的 时 


发 生 的 干旱 事件 数 (次 'a ) ;一 次 干旱 事件 的 历时 和 
强度 分 别 指 干旱 事件 持续 月 数 和 干旱 事件 中 SPEI 
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均值 的 绝对 值 , 本 研究 中 干旱 历时 和 干旱 强度 通过 对 
20 a 期 间 各 干旱 事件 的 历时 和 强度 取 均 值得 到 。 基 
于 干旱 特征 变量 构建 干旱 危险 性 指数 (DHI)” ,计算 
4 种 情景 下 各 模式 DHI, 再 对 4 种 情景 的 各 模式 结 
取 均 值 。DHI 计 算 公 式 如 下 : 
DHI = W(F x DY + W,(F x I)' (1) 
XB LOW. LWaZA BUS 28786, 15] 79 0.5 (FX D)' (Fx1I)' 分 
别 为 干旱 频次 (F) 与 干旱 历时 (D) FIE GO TH 
乘 所 得 乘积 的 归 一 化 值 , 归 一 化 处 理 使 用 极 值 标准 
化 法 ™。 
(2) 暴露 度 和 脆弱 性 分 析 
本 研究 选用 人 口 数量 .CDP 和 耕地 占 比 3 个 正 
问 指 标 对 烘 露 度 进行 量化 ,对 3 个 指标 进行 归 一 化 
处 理 后 取 等 权重 均值 ,得 到 干旱 暴露 度 指数 
(DEI)。 选 用 人 均 GDP 和 灌溉 占 比 2 个 负 癌 指标 对 
脆弱 性 进行 量化 ,对 2 个 指标 进行 归 一 化 处 理 后 取 
等 权重 均值 ,得 到 干旱 脆弱 性 指数 (DVI)。 
由 于 人 口 和 GDP 预 佑 数据 是 基于 2010 年 的 数 
据 生成 ,因此 对 于 人 口 数量 .CDP、 人均 CDP, 基 准 期 
指标 使 用 2010 年 数据 ,而 未 来 温 升 情景 下 指标 则 使 
用 20 a 数据 均值 。 对 于 耕 地 占 比 和 灌溉 占 比 2 个 指 
标 ,4 种 情景 下 指标 均 使 用 20 a 数据 均值 。 
(3) 旱灾 风险 指数 计算 
将 3 个 综合 指数 DHI、DEI、DVI 相 乘 ,得 到 4 种 
情景 的 旱灾 风险 指数 (DRI) ,计算 公式 如 下 : 
DRI = DHI x DEI x DVI (2) 
1.3.4 DRI 2 i4 46.8 5 ELAH 
(1) 空间 自 相 关 分 析 
基于 GeoDa 软件 对 4 种 情景 的 DRI 进 行 空 间 自 


日 变量 则 为 8 个 指标 (干旱 频次 干旱 历时 干旱 强 
BE A O Zt GDP 耕地 占 比 ,人均 CDP 灌溉 占 比 ) 
的 变化 , 即 每 个 温 升 情景 下 指标 值 与 基准 期 指标 值 
的 差 值 。 利 用 R 语 言 CD 包 实现 自 变 量 最 优 离散 化 
处 理 并 执行 地 理 探测 器 分 析 ”。 


2 结果 与 分 析 


2.1 气候 模式 模拟 性 能 评估 

由 于 SPEI 基 于 降水 量 和 PET 计 算得 到 ,CMIP6 
数据 对 中 国 降水 量 和 PET 模 拟 的 准确 性 可 反映 其 
对 干旱 的 模拟 能 力 。 利 用 1961 一 2014 年 的 CN05.1 
观测 数据 对 CMIP6 降水 和 PET 模 拟 数据 进行 了 评 
估 ,气候 模式 的 模拟 性 能 在 标准 化 泰勒 图 上 表现 为 
距离 观测 值 越 近 越 好 (图 2)。 降 水 量 模拟 值 与 观测 
值 的 + 范围 为 0.52~0.71,PET 的 + 范围 为 0.88~0.93。 
各 模式 降水 量 的 SD 均 大 于 1, 最 大 值 为 1.61, 而 PET 
的 SD 范围 为 0.93~1.39。 降 水 量 和 PET 的 RMSE 范 
围 分 别 为 0.78~1.32 和 0.37~0.68。 总 体 来 看 ,各 模 
式 对 降水 量 和 PET 的 模拟 性 能 良好 ,可 用 于 分 析 中 
国 干 旱 危 险 性 。 
2.2 温 升 情景 确定 结果 

统计 20 个 气候 模式 达到 2% 、 3 C4 MC 温 升水 
HY 20 a 时 间 段 ( 表 2) 发 现 ,所 有 气候 模式 都 能 达 
到 2 *C.3 % 温 升水 平 , 其 中 有 14 个 模式 达到 4 CT. 
升水 平 。 
23 DHI、DEI、DVI 3 个 综合 指数 计算 结果 

计算 4 种 情景 下 的 DHI、DEI、DVI, 使 用 分 位 数 
法 分 别 将 3 个 综合 指数 值 分 为 5 个 等 级 ,得 到 3 个 综 
合 指数 等 级 图 (图 3) ,第 1 等 级 到 第 5 等 级 (G1~G5) 


相关 性 分 析 。 采 用 全 局 Moran' s 1 从 整个 研究 区 尺 
度 检测 空间 自 相关 性 ,指数 的 取 值 范 围 为 -1~1, 大 
于 0 表示 正 相 关 , 小 于 0 表示 负 相 关 , 等 于 0 则 表示 
随机 分 布 。 采 用 局 部 Moran’ s 1 作为 局 部 空间 关联 
HERCLISA) ,分 析 空 间 集聚 特征 ”。 

(2) 地 理 探 测 需 

地 理 探 测 顺 是 探测 空间 分 异性 及 其 驱动 因素 
的 工具 ,包括 因子 探测 .交互 作用 探测 .风险 区 探测 
和 生态 探测 呈 。 本 文 基于 地 理 探测 需 分 析 旱 灾 风 
险 变 化 空间 分 异性 的 影响 因子 ,使 用 因子 探测 和 交 
互 作用 探测 得 到 衡量 自 变 量 对 因 变 量 解释 度 的 4 
值 。 将 每 个 温 升 情景 下 DRI 与 基准 期 DRI 作 差 , 得 
到 旱灾 风险 变化 ,作为 输入 地 理 探 测 需 的 因 变 量 。 


所 对 应 的 指数 值 依 次 增 大 。DHI 等 级 呈现 西北 与 东 
南 相 对 较 高 的 格局 ,与 中 国 西北 干旱 历时 长 .强度 
大 以 及 中 国 东南 干旱 频次 高 的 特点 相对 应 。DHI 
整体 上 随 着 温 升 水 平 的 增加 而 增 大 ,表明 未 来 中 国 
整体 呈现 暧 干 化 趋势 。DEI 等 级 呈 东 高 西 低 的 格局 
是 整体 随时 间 的 推移 而 增加 , DVI ERN E 
低 的 格局 且 整 体 随时 间 的 推移 而 降低 ,这 与 中 国 社 
会 经 济 发 展 水 平 东 高 西 低 的 格局 相符 ,未 来 由 于 人 
口 增长 和 经 济 发 展 , 暴 露 度 有 所 提高 而 脆弱 性 有 所 
降低 。 
2.4 DRI 空 间 格局 分 析 

使 用 分 位 数 法 将 DRI 数 值 分 为 5 个 等 级 ， 
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(a) 降水 量 
l 0.00 959 


观测 什 观测 值 
十 ACCESS-CM2  ® ACCESS-ESMI-5 x CanESMS E CanESM5-1 A CAS-ESM2-0 
© CMCC-ESM2 米 EC-Earth3 十 EC-Earth3-CC © EC-Earth3-Veg X EC-Earth3-Veg-LR 
FGOALS-g3 A FIO-ESM-2-0 Ô GFDL-ESM4 米 INM-CM4-8 十 INM-CM5-0 


© IPSL-CM6A-LR X MIROC6 © MPIESMI-2-HR A MPLIESMI2-LR ©@ MRI-ESM2-0 


IÈ :r, SD, RMSE PET 分 别 为 相关 系数 ,标准 差 . 均 方 根 误差 潜在 蒸 散 量 。 
图 2 各 气候 模式 降水 量 和 PET 的 标准 化 泰勒 图 


Fig. 2 Standardized Taylor diagrams of precipitation and PET for each climate model 


表 2 各 气候 模式 2%C.3 %C.4s%C 温 升 情景 年 份 


Tab.2 Years for 2 ^C, 3 C, and 4 °C temperature rise scenarios for each climate model 


"m 年 从 
TERR 2 CHIH 3 CIIM RE 4'CIHJHIR 

ACCESS-CM2 2024—2043 2043—2062 2058—2077 
ACCESS-ESMI-5 2029—2048 2050—2069 2068—2087 
CanESM5 2021—2040 2037—2056 2050—2069 
CanESM5-1 2021—2040 2036—2055 2050—2069 
CAS-ESM2-0 2029—2048 2047—2066 2063—2082 
CMCC-ESM2 2033—2052 2049—2068 2064—2083 
EC-Earth3 2026—2045 2048—2067 2064—2083 
EC-Earth3-CC 2039—2058 2056—2075 2069—2088 
EC-Earth3-Veg 2031—2050 2049—2068 2065—2084 
EC-Earth3-Veg-LR 2034—2053 2054—2073 2067—2086 
FGOALS-g3 2039—2058 2066—2085 一 

FIO-ESM-2-0 2031—2050 2050—2069 2066—2085 
GFDL-ESMA 2037—2056 2061—2080 一 

INM-CM4-8 2037—2056 2060—2079 一 

INM-CM5-0 2035—2054 2063—2082 - 

IPSL-CM6A-LR 2029—2048 2045—2064 2060—2079 
MIROC6 2037—2056 2062—2081 2080—2099 
MPI-ESMI-2-HR 2039—2058 2064—2083 一 

MPI-ESMI-2-LR 2040—2059 2063—2082 一 

MRI-ESM2-0 2025—2044 2052—2071 2071—2090 


注 : -表示 气候 模式 不 能 达到 4 "CIRCE o 
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(h) 3 °C 温 升 情景 DEI 
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(f) 2 °C 温 升 情景 DVI 


G) 3 “C 温 升 情景 DVI 
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注 :DHI、DEI、DVI 分 别 为 干旱 危险 性 指数 .干旱 暴露 度 指数 .干旱 脆弱 性 指数 。 下 同 。 
图 3 基准 期 和 2% .3% .4% 温 升 情 景 下 中 国 DH DEI. DVI 等 级 分 布 
Fig.3 Distributions of DHI, DEI, and DVI grades for the baseline period and the 2 °C, 3 °C, and 4 °C 


temperature rise scenarios in China 


Gl~G5 所 对 应 的 DRI 数 值 依次 增 大 (图 4a.d\g、j)。 
DRI 等 级 呈 东 高 西 低 的 格局 , 且 随 温 升水 平 的 增加 
而 整体 增 大 。 全 局 空间 自 相关 分 析 结 果 显 示 ,4 种 
情景 下 的 全 局 Moran’ s 1 值 分 别 为 0.851、0.864、 
0.860 和 0.855, 且 都 通过 了 显著 性 检验 (P<0.001)。 
根据 Moran 散 点 图 (图 4c fi 1) ,大 多 数 散 点 均 落 在 
第 一 象限 (高 高 区 ) 和 第 三 象限 ( 低 低 区 ) ,少数 散 点 
落 在 了 第 二 象限 ( 低 高 区 ) 和 第 四 象限 (高 低 区 ) , 表 
HH DRI 空 间 分 布 呈 现 以 高 值 集聚 和 低 值 集聚 为 主 的 


空间 正 相 关 。 进 一 步 绘制 了 0.05 显著 性 水 平 下 的 
LISA 集聚 图 (图 4b ,e .h Kk) ,将 所 有 格 点 分 为 5 种 类 
型 :不 显著 高 高 区 ARREK 、 低 高 区 和 高 低 区 。 高 
值 集聚 区 主要 分 布 在 地 区 A 、B、C.D, 而 低 值 集聚 区 
主要 分 布 在 地 区 EF、G。 

全 国 及 7 个 自然 地 区 的 DRI 箱 线 图 (图 5) 显 示 : 
从 全 国 来 看 ,DRI 的 最 大 值 .中 位 数 和 均值 从 基准 期 
到 4% 温 升 情景 依次 增 大 ,3 个 温 升 情景 下 的 极 差 依 
次 增 大 ,平均 值 全 部 大 于 中 位 数 。 分 区 域 来 看 ,在 
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注 :DRI 为 旱灾 风险 指数 ;LISA 为 局 部 空间 关联 指标 。 下 同 。 
图 4 DRI 等 级 分 布 及 空间 自 相 关 性 分 析 结 


Fig. 4 DRI grade distributions and results of spatial autocorrelation analysis 


基准 期 ,DRI 的 最 大 值 出 现在 C 地 区 ,中 位 数 与 均值 
的 排序 一 致 ,各 地 区 由 高 到 低 依次 为 :B、C、A、D、E、 
FG; 在 3 个 温 升 情景 下 ,DRI 最 大 值 仍 出 现在 C 地 
区 ,中 位 数 与 均值 最 高 的 3 个 地 区 为 地 区 BC.D, 最 
低 的 地 区 为 G; 整 体 来 看 ,7 个 自然 地 区 DRI 随 着 温 
升水 平 的 增加 而 增 大 。 
2.5 旱灾 风险 变化 分 析 

将 2% .3%,、4% 温 升 情景 下 的 旱灾 风险 变化 


值 划分 为 6 个 等 级 ,得 到 空间 分 布 如 图 6 所 示 。 整 
体 来 看 ,3 个 温 升 情景 下 的 DRI 相 对 于 基准 期 以 增 
大 为 主 ,在 中 国 东部 沿海 地 区 增 大 较 明 显 。 随 着 温 
升水 平 的 升 高 ,DRI 增 加 值 整体 增 大 。 

因子 探测 和 交互 作用 探测 结果 如 图 7 所 示 。 就 
单 因 子 而 言 ,在 3 个 温 升 情景 下 ,8 个 影响 因子 的 g 
值 均 通过 了 0.001 水平 的 显著 性 检验 ,g 值 排名 前 3 
的 影响 因子 均 为 人 口 数量 变化 .CDP 变 化 .耕地 占 
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图 5 全 国 和 7 个 自然 地 区 DRI 箱 线 图 
Fig. 5 Box plots of DRI for the whole country and 7 natural regions 
(a) 2 "C 温 升 情景 (b) 3 "C 温 升 情景 (c) 4 °C 温 升 情景 
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Fig. 6 Distributions of changes in drought risk under 2 *C, 3 °C, and 4 °C temperature rise scenarios 


(a) 2 "C 温 升 情景 (b) 3 "C 温 升 情景 (c) 4 "C 温 升 情景 


G XXX X X X% X, Xs X; X; Xe X X X, X, X, X, X, X, 
影响 因子 影响 因子 影响 因子 


IE FE 
qfi 0.0 01 02 03 04 05 06 0.7 0.8 09 


W:X A A Aa X AX AX ADA FFR FRA PRE A O ži GDP ,耕地 占 比 人均 CGDP 灌溉 占 比 的 变化 ; 
4 为 自 变 量 对 因 变 量 解释 度 。 
图 7 因子 探测 和 交互 作用 探测 结果 


Fig.7 Results of factor detection and interaction detection 
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比 变化 ;在 2 % 温 升 情 景 下 ,耕地 占 比 变化 对 旱灾 风 
险 变 化 的 解释 度 最 大 ,为 61%; 在 3 % .4% 温 升 情 景 
下 ,GDP 变 化 对 旱灾 风险 变化 的 解释 度 最 大 ,为 
63%。 就 影响 因子 交互 作用 而 言 ,交互 作用 类 型 为 
双 因 子 增强 和 非 线 性 增强 ;在 3 个 温 升 情景 下 ,对 早 
灾 风 险 变 化 解释 力 最 大 的 2 对 交互 作用 影响 因子 均 
为 耕地 占 比 变化 人口 数量 变化 .耕地 占 比 变化 门 
GDP 变 化 ,其 交互 作用 对 风险 变化 的 解释 程度 在 
75% 以 上 ,表明 承 灾 体 暴 露 度 变化 对 旱灾 风险 变化 
的 解释 能 力 较 强 。 综 合 来 看 ,影响 因子 两 两 交互 作 
用 相对 于 单 因子 而 言 能 更 好 地 解释 旱灾 风险 变 
化 。 王 旱 频 次 变化 .干旱 历时 变化 .干旱 强度 变化 
的 单 因 子 g 值 及 两 两 交互 作用 的 g 值 在 2 C3 C, 
4 人 温 升 情景 下 依次 增 大 ,表明 干旱 危险 性 变化 对 
旱灾 风险 变化 的 影响 随 着 温 升 加 剧 而 更 加 明显 。 


3 讨论 


目前 ,基于 未 来 情景 数据 进行 中 国 旱 灾 风 险 预 估 
的 研究 相对 较 少 , 且 已 有 人 研究 多 为 区 域 性 研究 。 
本 文 引 入 未 来 预 估 数据 所 得 中 国 旱灾 风险 的 空间 
分 布 与 赵 佳 琪 等 ”基于 历史 数据 所 得 结果 类 似 , 表 
明了 研究 结果 的 可 乱 性 。 与 前 人 研究 相 比 ,本 文 从 
不 同 温 升 情景 的 角度 开展 旱灾 风险 评估 ,并 分 析 了 
未 来 旱灾 风险 变化 的 驱动 因素 ,有 助 于 更 好 地 认识 
中 国 旱灾 风险 随 着 全 球 温 升 加 剧 而 发 生 的 变化 ,为 
气候 变化 风险 防范 提供 了 参考 依据 。 然 而 ,本 研究 
还 存在 一 些 不 足 之 人 处。 由 于 未 来 气候 变化 和 人 类 
活动 的 复杂 性 以 及 气候 模式 的 不 确定 性 ,基于 
SSP5-8.5 情 景 的 气候 模式 数据 所 得 结果 具有 较 大 
的 不 确定 性 ,中国 未 来 气候 是 否 持续 暖 干 化 还 有 和 
进一步 探讨 。 同 时 ,指标 选取 数据 处 理 方式 .权重 
设置 等 也 会 造成 旱灾 风险 评估 结果 的 不 确定 性 。 
此 外 ,本 文 对 中 国 所 有 格 点 进行 了 旱灾 风险 变化 的 
全 局 驱动 力 探测 ,但 不 同 地 区 由 于 自然 和 社会 经 济 
条 件 存在 差异 ,旱灾 风险 变化 的 主导 因素 可 能 有 所 
不 同 , 未 来 可 进一步 探测 不 同 地 区 的 局 域 驱动 力 并 
分 析 其 差异 。 


4 结论 


本 研究 利用 20 个 CMIP6 气 候 模 式 数 据 和 承 灾 
体 预 佑 数据 ,对 致 灾 因子 危险 性 、 承 灾 体 暴 露 度 和 


脆弱 性 进行 量化 ,计算 了 DHI、DEI 和 DVI, 从 而 将 3 
个 综合 指数 相 乘 得 到 DRI, 评 估 了 2 % 、 3 C4 "Cil 
升 情景 下 的 中 国 旱灾 风险 及 其 相对 于 基准 期 的 变 
化 ,主要 结论 如 下 : 

(1) DHI, DEI FU DVI Æ £ BH 5 B5) zs [8] 4) 4g 
律 。DHI 整 体 上 表现 为 西北 和 东南 相对 较 高 ,DFI 
和 DVI 则 分 别 表现 为 东 高 西 低 西高 东 低 。 随 着 温 
升水 平 的 升 高 ,干旱 危险 性 增加 明显 , 承 灾 体 暴露 
度 整 体 有 所 增加 ,而 承 灾 体 脆弱 性 有 所 降低 。 

(2) DRI 呈 现 明 显 的 空间 集聚 特征 ,高 值 空间 
集聚 区 主要 分 布 在 东部 季风 区 , 低 值 空间 集聚 区 主 
要 分 布 在 西部 地 区 ,中 国 旱 灾 风 险 在 基准 期 和 不 同 
温 升 情景 下 都 呈 东 高 西 低 的 特点 。 

(3) 与 基准 期 相 比 ,3 个 温 升 情景 下 的 DRI 以 增 
大 为 主 ,中 国 东部 沿海 地 区 增 大 最 为 明显 。 随 着 温 
升水 平 的 升 高 ,风险 增加 值 整体 呈 增 长 趋势 。 通 过 
对 旱灾 风险 变化 进行 全 局 驱动 力 探测 ,发现 人 口 数 
量变 化 .GDP 变化 和 耕地 占 比 变化 是 主要 驱动 因 
T ,因子 两 两 交互 作用 相对 于 单 因子 对 旱灾 风险 变 
化 有 更 强 的 解释 力 。 
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Abstract: Drought is one of the most disastrous extreme climate events. Studying the changes in drought risk 
against the background of future global warming is beneficial for scientifically advancing disaster prevention and 
reduction work deployment. The standardized precipitation evapotranspiration index was calculated using data 
from 20 climate models from the sixth phase of the Coupled Model Intercomparison Project. Drought characteris- 
tic variables were extracted for the baseline period and global temperature rise scenarios of 2 °C, 3 °C, and 4 *C in 
China, and the drought hazard index (DHI) was calculated. Based on the disaster- bearing body projection data, 
the drought exposure index (DEI), drought vulnerability index (DVI), and drought risk index (DRI) were compre- 
hensively calculated. The distribution pattern of drought risk in China was further analyzed, and a spatial attribu- 
tion analysis of future drought risk changes was performed using a geodetector. The results showed that the spa- 
tial distribution of DHI, DEI, and DVI exhibited higher values in the northwest and southeast, a pattern of being 
high in the east and low in the west, and a trend of being high in the west and low in the east, respectively. Based 
on this, the DRI specially showed a spatial positive correlation dominated by high- and low- value clustering. 
With an increase in the temperature rise level, the future drought risk will mainly increase across China, and the 
increase in the eastern coastal areas would be the most obvious. Changes in population, gross domestic product, 
and the proportion of cultivated land were found to be the main factors affecting changes in drought risk. 


Key words: drought risk; CMIP6; temperature rise scenario; spatial autocorrelation; geodetector 


